Decay constants and radiative decays of heavy mesons in light-front quark model 
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We investigate the magnetic dipole decays V — * P7 of various heavy-flavored mesons such as 
(D, D* , D 3 , D* ,r/ c , J/tp) and (B, B* , B a , B* ,rjb, T) using the light-front quark model constrained by 
the variational principle for the QCD-motivated effective Hamiltonian. The momentum dependent 
form factors Fvp(q 2 ) for V — » P7* decays are obtained in the q + = frame and then analytically 
continued to the timelike region by changing q x to zqx in the form factors. The coupling constant 
gvp-y for real photon case is then obtained in the limit as q 2 — > 0, i.e. gvp-y = Fyp (q 2 = 0). The 
weak decay constants of heavy pseudoscalar and vector mesons are also calculated. Our numerical 
results for the decay constants and radiative decay widths for the heavy-flavored mesons are overall in 
good agreement with the available experimental data as well as other theoretical model calculations. 



I. INTRODUCTION 

The physics of exclusive heavy meson decays has pro- 
vided very useful testing ground for the precise deter- 
mination of the fundamental parameters of the standard 
model(SM) and the development of a better understand- 
ing of the QCD dynamics. While the experimental tests 
of exclusive heavy meson decays are much easier than 
those of inclusive one, the theoretical understanding of 
exclusive decays is complicated mainly due to the nonper- 
turbative hadronic matrix elements entered in the long 
distance nonperturbative contributions. Since a rigrous 
field-theoretic formulation with a first principle applica- 
tion of QCD to make a reliable estimates of the nonper- 
turbative hadronic matrix elements has not so far been 
possible, most of theoretical efforts have been devoted to 
looking for phenomenological approaches to nonpertur- 
bative QCD dynamics. 

Along with various exclusive processes such as leptonic, 
semileptonic and rare decays of heavy mesons, the one- 
photon radiative decays from the low-lying heavy vec- 
tor(V) to heavy pseudoscalar(P) mesons, i.e. magnetic 
dipole U(l 3 Si) — > P(l 1 5o)7 transitions, have been con- 
sidered as a valuable testing ground to pin down the 
best phenomenological model of hadrons. For example, 
the calculations of D* — > D-f and B* — » Bj radia- 
tive decays have been investigated by various theoreti- 
cal approaches, such as the quark model [l|, [H, H, 0, @, 
light cone QCD sum rules H, 0], heavy quark effec- 
tive theory (HQET) cloudy bag model(CBM) QjJ, 
and chiral perturbation theory [Til ]. Recently, the ra- 
diative decays between two heavy quarkonia such as 
J/ip — ► r] c j and T — > 77^,7 have also been studied by 
the potential nonrelativistic QCD(pNRQCD) [l2| and 
relativistic quark model (l3l [lij approaches. In our 
previous light-front quark model(LFQM) analysis [HI, 
WL [ItJ based on the QCD-motivated effective Hamil- 
tonian, we have analyzed various exlcusive processes 
such as the 0~ — > 0~ semileptonic heavy meson de- 
cays [l6|], rare B — > Kl + l~{l = e,/x,r) decays [l7| 
and radiative V — > P7 and P — * U7 decays of light- 
flavored mesons(7r, p, ui, K, K* , </>, 77, 7/) [l5|] and found a 
good agreement with the experimental data. 



The main purpose of this work is to investigate the 
magnetic dipole transition V — > P7 for the heavy- 
flavored mesons such as (D, D* , D s ,D * r] C: J_H >) and 
(B,B*,B s ,B*,ri b ,T) using our LFQM [Il,[lE[i3]- Since 
the experimental data available in this heavy-flavored 
sector are scanty, predictions of a model, if found reliable, 
can be utilized quite fruitfully. In addition, we calculate 
the weak decay constants of heavy pseudoscalar and vec- 
tor mesons, which play important roles in many aspects, 
such as in the determination of Cabibbo-Kobayashi- 
Maskawa matrix elements, in the leptonic or nonleptonic 
weak decays of mesons, and in the neutral D — D or 
B - B mixing process, etc. Our LFQM [H, HE (13 use d 
in the present analysis has a couple of salient features 
compared to other LFQM 0, [l3j] analysis: (1) We have 
implemented the variational principle to QCD-motivated 
effective LF Hamiltonian to enable us to analyze the me- 
son mass spectra and to find optimized model parame- 
ters, which are to be used subsequently in the present 
investigation. Such an approach can better constrain 
the phenomelogical parameters and establish the extent 
of applicability of our LFQM to wider ranging hadronic 
phenomena. (2) We have performed the analytic contin- 
uation from the spacelike(<7 2 < 0) region to the physical 
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- M P ) 2 ) to obtain the 
V — * P7* transitions. 



timelike region(0 < q 
decay form factors Fyp 
The Drell-Yan-West(<7 + = q° + q 3 = 0) frame is useful 
because only valence contributions are needed, i.e. the 
hadronic matrix element (P| J^lV) is represented as the 
overlap of valence wave function, as far as the " fi = +" 
component of the current is used. 

The paper is organized as follows: In Sec. II, we briefly 
describe the formulation of our LFQM [H Hi] and the 
procedure of fixing the model parameters using the vari- 
ational principle for the QCD-motivated effective Hamil- 
tonian. The decay constants and radiative V — > P7 de- 
cay widths for heavy-flavored mesons are then uniquely 
determined in our model calculation. In Sec. Ill, the for- 
mulae for the decay constants of pseudoscalar and vector 
mesons as well as the decay widths for V — > P7 in our 
LFQM are given. To obtain the g 2 -dependent transition 
form factors Fyp(q 2 ) for V — > P7* transitions, we use 
the q + = frame(i.e. q 2 = — < 0) and then analyti- 
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cally continue the spacelike results to the timelike q 2 > 
region by changing to zqj_ in the form factor. The 
coupling constants gvp-y needed for the calculations of 
the decay widths for V — > P"f can then be determined in 
the limit as q 2 — ► 0, i.e. gypj — Fyp(q 2 — 0). In Sec. 
IV, we present our numerical results and compare with 
the available experimental data as well as other theoret- 
ical model predictions. Summary and conclusions follow 
in Sec.V. 



The covariant forms of the spin-orbit wave functions 
for pseudoscalar and vector mesons are given by 
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-u(pi,Ai) 
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(4) 



II. MODEL DESCRIPTION 



where Mq = \/M^ — (toi — m^) 2 and M§ is the invariant 
meson mass square Mq defined as 



The key idea in our LFQM 0, [3 

for mesons is 

to treat the radial wave function as trial function for 
the variational principle to the QCD-motivated effective 
Hamiltonian saturating the Fock state expansion by the 
constituent quark and antiquark. The QCD-motivated 
Hamiltonian for a description of the ground state meson 
mass spectra is given by 



JJz 



to| + k 2 



I nlrni 



= [H + V q g]\* 



' ntm/ 



(1) 



where k = (k j_ ,k z ) is the three-momentum of the con- 
stituent quark, M qq is the mass of the meson, and \^ n i^) 
is the meson wave function. In this work, we use two in- 
teraction potentials V qq for the pseudoscalar(0 '") and 
vector(l ) mesons: (1) Coulomb plus harmonic osclla- 
tor(HO), and (2) Coulomb plus linear confining poten- 
tials. In addition, the hyperfine interaction, which is es- 
sential to distinguish vector from pseudoscalar mesons, 
is included for both cases, viz., 



^k 2 

^o 2 = E— 



rrii 



(5) 



The polarization vectors e^{J z ) of the vector meson with 
four momentum P are given by 



e*(±l) 



0,— e x (±)-Px,ex(±l) 



V2 ' 



P 4 



p2 _ M 2 

p+ ' J 



(6) 



The spin-orbit wave functions satisfy the following rela- 
tions 



^ ; "^"Ai A2 ^Ai A2 ^' 



(7) 



Ai A2 



for both pseudoscalar and vector mesons. For the radial 
wave function </>, we use the same Gaussian wave function 
for both pseudoscalar and vector mesons 



Aq, 2 S S 

V qq = Vo+Vhyp = a+Vconf-^ + r^— -^Wcoui, (2) 

3r 3 m q m q 

where V CO nf = br(r 2 ) for the linear(HO) potential and 
(Sq-Sg) = l/4(— 3/4) for the vector (pseudoscalar) meson. 

The momentum space light-front wave function of the 
ground state pseudoscalar and vector mesons is given by 



*i r 5&(aii > ku.,Ai) =^ 1 J J a (a;i > k 4 x)^(a; i ,kix), (3) 



where (f>(xi,\s.i±) is the radial wave function and ^ 
is the spin-orbit wave function, which is obtained by the 
interaction independent Melosh transformation from the 
ordinary equal-time static spin-orbit wave function as- 
signed by the quantum numbers J PC . The model wave 
function in Eq. ([3]) is represented by the Lorentz- invariant 
variables, Xi — p\ / ' P + , kij_ = Pi± — XiP± and A$, where 
pf > = (p+ ) p- ) p x ) = (P° + P 3 ,(A/ 2 + P 2 L )/P+,P ± ) 
is the momentum of the meson M, p^ and A^ are the 
momenta and the helicities of constituent quarks, respec- 
tively. 



4^3/4 I ah. 

cj>(x i ,k iX ) = ^j T ^/^exp(-k 2 /2p% (8) 

where (3 is the variational parameter. When the longi- 
tudinal component k z is defined by k z = (x — 1/2) Mq + 
(to 2 — mf)/2Mo, the Jacobian of the variable transfor- 
mation {x, kx} —tk= (kj_, k z ) is given by 



dx 



M 



1 



Ml 



(9) 



4X1X2 

Note that the free kinetic part of the Hamiltonian Hq = 



2 + k 2 + y m| + k 2 is equal to the free mass operator 
Mo in the light-front formalism. 

The normalization factor in Eq. ([8|) is obtained from 
the following normalization of the total wave function, 



'dx [ ^\^{x,ku.)\ ! 



1. 



(10) 



Our variational principle to the QCD-motivated effec- 
tive Hamiltonian first evaluate the expectation value of 
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TABLE I: The constituent quark mass[GeV] and the Gaussian 
paramters /3[GeV] for the linear and HO potentials obtained 
by the variational principle, q = u and d. 



Model m q 


m 3 


m c 


m b 




/3sc 


/3cc 


Pqb 


Psb 


Pbb 


Linear 0.22 


0.45 


1.8 


5.2 


0.468 


0.502 


0.651 


0.527 


0.571 


1.145 


HO 0.25 


0.48 


1.8 


5.2 


0.422 


0.469 


0.700 


0.496 


0.574 


1.803 
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n t .(29SO 
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(470) 


K(494) 


(478) 




0X7821 






71(140) 




K*(892) 
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(Exp.) 
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FIG. 1: (color online). Fit of the ground state meson 
masses [Me V] with the parameters given in Table[I] The (p, tt), 
(77,7/), and (ijJ,(f>) masses are our input data. The masses of 
(oj — <j>) and (77 — 77') were used to determine the mixing angles 
of W — <j> and 77 — 77' [l5j|. respectively. 



the central Hamiltonian iJo + Vo, i.e. (<^|(i?o + Vb)|0) 
with a trial function k^j_) that depends on the vari- 
ational parameters (3 and varies (3 until (0|(i?o + Vb) |0) 
is a minimum. Once these model parameters are fixed, 
then, the mass eigenvalue of each meson is obtained by 
M q q — {4>\{Hq + V q q)\4>) . On minimizing energies with re- 
spect to (3 and searching for a fit to the observed ground 
state meson spectra, our central potential Vo obtained 
from our optimized potential parameters (0 = —0.72 
GeV, b = 0.18 GeV 2 , and a s = 0.31) [H for Coulomb 
plus linear potential were found to be quite comparable 
with the quark potential model suggested by Scora and 
Isgur [Hi where they obtained a = -0.81 GeV, b = 0.18 
GeV 2 , and a s = 0.3 ~ 0.6 for the Coulomb plus linear 
confining potential. More detailed procedure of deter- 
mining the model parameters of light and heavy quark 
sectors can be found in our previous works [Hill. Our 
model parameters (m, /?) for the heavy quark sector ob- 
tained from the linear and HO potential models are sum- 
marized in Table [U 

Our predictions of the ground state pseudoscalar and 



vector meson mass spectra obtained from the linear po- 
tential parameters were already shown in [l6l |. In this 
work, we include the results obtaintcd from the HO po- 
tential parameters as well and summarize them in Fig. [1] 
Our predictions of the ground state meson mass spectra 
obtained from both linear and HO parameters agree with 
the experimental data [l!| within 6% error. We should 
note that our previously predicted mass of B*, Mb* = 
5424[5471] MeV obtained from the linear[HO] parameters 
is in good agreement with the very recent CLEO data, 
M B > = 5414±l±4MeV[H andM B « = 5411.7±1.6±0.6 
MeV [2l| • For the experimentally unmeasured mass of 1% 
meson, our prediction of mass difference between the two 
bottomonia Am(= M r - M m ) = 34[263] MeV obtained 
from the linear [HO] parameters is consistent with cur- 
rent theoretical estimates (from pcrt urbative QCD and 
lattice NRQCD), Am = 34 ~ 141 MeV Q. As we shall 
see in our numerical calculations, the radiative decay of 
T — * r]bl might be useful to determine the mass of 77;, 
experimentally since the decay width T(T — * 77^7) is very 
sensitive to the value of Am, viz. Y cx (Am) 3 . 



III. DECAY CONSTANTS AND RADIATIVE 
DECAY WIDTHS 

The decay constants of pseudoscalar and vector mesons 
are defined by 



(0\q^ l5 q\P) = if P P», 
{0\rfq\V(P,h)) = fvM v e^h) 7 



(11) 



where the experimental value of vector meson decay con- 
stant fy is extracted from the longitudinal(/i = 0) polar- 
ization. In the above definitions for the decay constants, 
the experimental values of pion and rho meson decay 
constants are / w w 131 MeV from it — > [w and f p w 220 
MeV from p — > e + e~ . 

Using the plus component (/i = +) of the current, one 
can easily calculate the decay constants. The explicit 
forms of pseudoscalar and vector meson decay constants 
are given by 



fp 
fv 



2a/6 
2V6 



dx d 2 k± 



A 



16/T 3 

dx d 2 k\ 



4>( x i k 



16/T 3 



ki 



(/>(x,k±), 
A - k 



Mo 



,(12) 



where A = xim\ + X\mi and A4q = Mq + mi + 777,2. 

In our LFQM calculation of V —> P-f decay process, we 
shall first analyze the virtual photon(7*) decay process 
so that we calculate the momentum dependent transi- 
tion form factor, Fvp(q 2 )- The lowest-order Feynman 
diagram for V — * P7* process is shown in Fig. [2] where 
the decay from vector meson to pseudoscalar meson and 
virtual photon state is mediated by a quark loop with 
flavors of consituent mass mi and 777,2. 
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FIG. 2: Lowes-order graph for V — > P7* transitions. 



The transition form factor Fyp(q 2 ) for the magnetic 
dipole decay of vector meson V(P) — ► P(P')-f*(q) is de- 
fined as 

(P(P')I J£ n | W ft)) = iee^ CT e,(P, ^P.Pyp^ 2 ), 

(13) 

where the antisymmetric tensor e^ pa assures electro- 
magnetic gauge invariance, q = P — P' is the four mo- 
mentum of the virtual photon, e„(P, h) is the polarization 
vector of the initial meson with four momentum P and 
helicity h. The kinematically allowed momentum trans- 
fer squared q 2 ranges from to g^ ax = (My — Mp) 2 . 

The decay form factor Fvp(q 2 ) can be obtained in the 
q + = frame with the "good" component of currents, 
i.e. (/, = +, without encountering zero- mode contribu- 
tions [23j. Thus, we shall perform our LFQM calculation 
in the q + = frame, where q 2 = q + q~ — = — qj_ < 0, 
and then analytically continue the form factor Pyp(q^) 
in the spacelike region to the timelike q 2 > region by 
changing q^ to iq± in the form factor. 

The quark momentum variables for V(q±q~2) — > P((fiQ2) 
transitions in the q + = frame are given by 



Pi 

PlX 



,T 1 P+, i 3+ = x 2 P+, 

XiPj. + kj_, p 2 J_ = X 2 Px - kx, 
Xl P + ,p'+ = X 2 P + , 

x 1 P'x + k'x,p' 2 x = ^P'x-k'x, (14) 



where xi = x and X2 = 1 — x and the spectator quark re- 
quires that P2 = p' 2 + an d P2_l = p'ax- ^ n ^ ne calculations 
of the decay form factor Fyp{q 2 ), we use '+ '-component 
of currents and the transverse(/i = ±1) polarization. For 
the longitudinal(/i = 0) polarization, it is hard to extract 
the form factor since both sides of Eq. (|13[) are vanishing 
for any q 2 value. 

The hadronic matrix element of the plus current, 
(,/+) = (P(P')\J+JV(P,h = +)) in Eq. (33) is then 
obtained by the convolution formula of the initial and 
final state light-front wave functions: 



(J + ) 



f dx f 
^2 ee j J J rf 2 kj_^(x,k'_L)(/)(x,kj 



OQt^A'CPi) +«a(Px)_i 



AA 



TZ n - 
XX' 



(15) 



where k'x = k^— x 2 qj_ and ee 3 is the electrical charge for 
j-th quark flavor. Comparing with the right-hand-side 
of Eq. (HI, i.e. eP+ F VP {Q 2 )q R / V2 where q R = q x + 
iq y , we could extract the one-loop integral, I (mi, rri2, q 2 ), 
which is given by 



7(TOi,m 2 ,g 2 ) = J J d 2 k 



(x,k'j_)</>(x,k ± ) 
xiM Q M' 
(k ± -q ± ) 2 - 



(16) 



where the primed factors are the functions of final state 
momenta, e.g. Mq — Mq (x,kj_). 

Then, the decay form factor Fyp(q 2 ) is obtained as 

F V p(q 2 ) = ei/(mi,TO 2 ,<7 2 ) + e 2 I{m 2l m 1 ,q 2 ). (17) 

The coupling constant gvPj I0r real photon(7) case can 
then be determined in the limit as q 2 — > 0, i.e. <?yp 7 = 
Fyp(q 2 — 0). The decay width for V — > P7 is given by 



nr- f<?y P7 fc 7 , 



(18) 



where a is the fine-structure constant and fc 7 = (My — 
Mp)/2My is the kinematically allowed energy of the out- 
going photon. 



IV. NUMERICAL RESULTS 

In our numerical calculations, we use two sets of model 
parameters (m, 0) for the linear and HO confining poten- 
tials given in Table U to perform, in a way, a parameter- 
free-calculation of decay constants and decay rates for 
heavy pseudoscalar and vector mesons. Although our 
predictions of ground state heavy meson masses are over- 
all in good agreement with the experimental values, we 
use the experimental meson masses in the computations 
of the radiative decay widths to reduce possible theoreti- 
cal uncertainties. But in the case of 775 , for which experi- 
mental data is not yet available, we use the model mass as 
M Vb = 9353 ± 50 MeV, i.e. we use slightly broader range 
Am = My 
in 0. 



M m = 60 - 160 MeV than that reported 



In Table HI! we present our predictions for the charmed 
meson decay constants (fp, fp* ,f D . , fp* , f Ve , fj /,;,) to- 
gether with lattice QCD [H SU, QCD sum rules [H, 
relativistic Bethe-Salpeter (BS ) model [27|, rela- 
tivized quark model [281 ] . and other relativistic quark 
model(RQM) [2^] predictions as well as the available 
experimental data [lj| [3(| HH H3] • Note that we extract 
the experimental value (fj/^)exp = (416 ± 6) MeV from 
the data r oxp (J/^_-> e+e") = 5.55±0.14±0.02 keV [H 
and the formula 13311 



T(V 



e ' e 



47T a 2 



fycv, 



(19) 
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TABLE II: Charmed meson decay constants(in unit of MeV) obtained from the linear[HO] parameters. 





So 


Id* 


Sn s 




he 


Sj/i> 


Linear [HO] 


211(194] 


254 [228] 


248 [233] 


290 [268] 


326 [354] 


360(395] 




Oil J_ 1 A^~' 2 
Zll It l i ±_22 


Z4D zt ^U_2 


OQ1 _|_ 1 O+S 
ZOl zt 


079 -1- 1 + ^ 
Z ( Z Zt 1D_20 






QCD [25] 


201±3±17 




249 ± 3 ± 16 








Sum-rules [26] 


204 ± 20 




235 ± 24 








BS [27] 


230 ± 25 


340 ± 23 


248 ± 27 


375 ± 24 


292 ± 25 


459 ± 28 


QM [28] 


240 ± 20 




290 ± 20 








RQM [29] 


234 


310 


268 


315 






Exp. 


222.6 ± 16.7^3 4 [30J 




282 ± 16 ± 7 [31] 




335 ± 75 [32] 


416 ±6 [19] 



where cy = 4/9 for V — J/i/j. Our predictions 
for the ratios fnjfn = 1.18[1.20] and f Vc /fj/y = 
0.91(0.90] obtained from the linear[HO] parameters are 
in good agreement with the available experimental data, 
C/W/dW = 1.27 ± 0.12 ±0.03(preliminary) [H| and 
CAje//.//y>)exp. = 0.81 ±0.19 [H,[32|, respectively. Our 
result for the ratio /d*//d» = 1.14(1.18] obtained from 
the linear [HO] parameters is also consistent with the 
quenched lattice result, 1.11(3) [24| and the BS one, 
1.10 ± 0.06 (27[. Overall, our results for the charmed 
meson decay constants are in good agreement with other 
theoretical model calculations as well as the experimental 
data. 

In Table HTT1 we show our results for the bottomed me- 
son decay constantsf/f?, /b* ■ £g s , fs* , Svb j /t ) together 
with lattice QCD [If, QCD sum rules [H, 

[36| . BS model [271 ]. relativized quark model (28|, and 
RQM [29[ predictions as well as the available experi- 
mental data [n| H3]- Note that we extract the exper- 
imental value (/x)exp = (715 ± 5) MeV from the data 
r cX p(T -> e+e") = 1.340 ± 0.018 keV [3 and Eq. pjl 
with cy = 1/9 for V = T. Our results for the ra- 
tios /b,// b - 1.24(1.32] and /b ; // b . = 1.23(1.32] ob- 
tained from the linear [HO] parameters are quite com- 
parable with the recent lattice results, 1.20(3)(1) (34| 
and 1.22(±|) [U for / B; // s » and 1.17(4)±£ UM for 
/b* //b* • For the T meson decay constant, our prediction 
fr = 529(893] MeV obtained from the linear [HO] param- 
eters slightly deviates from the extracted experimental 
value (fr)exp = (715 ± 5) MeV. Other model calcula- 
tions for f r such as 498 ±20 MeV from the BS model H3 
and 836 MeV from effective Lagrangian satisfying heavy- 
quark spin symmetry(HQSS) [38| also show some devi- 
ations from the experimental value. Our result for the 
ratio f Vb /fr — 0.96(0.91] obtained from the linear[HO] 
parameters is to be compared with the /fx ~ 1 in 
HQSS limit [38[. For these heavy bottomed meson decay 
constants, we observe an overall agreement between our 
results and other theoretical ones. 

In Table IIV1 we present our results of the coupling 
constants gyp 7 (in unit of GeV -1 ) for radiative V — ► P7 
decays together with other QM calculations [l], H, Q as 
well as the available experimental data. The experimen- 
tal values for {gj/^ Vc7 ) cxp = 0.57 ± 0.11 for J/ip — > r; c 7 
and (3d-+d+ 7 )ox P = -(0.50±0.12) for D*+ -> D+7 pro- 



cesses are extracted from the branching ratios Br (J/ip — ► 
%7)cx P = (1.3 ± 0.4)% and Br(£>*+ -> P+ 7 )c Xp = 
(1.6±0.4)% together with the full widths of r tot (J/</>) = 
93.4±2.1keVandr tot (L>*+) = 96±22 keV The op- 
posite sign of coupling constants for D* + and D* + decays 
compared to the charmonium J/ip decay indicates that 
the charmed quark contribution is largely destructive in 
the radiative decays of D* + and D* + mesons. Similarly, 
we see that the bottomed quark contribution is largely 
destructive in the radiative decay of B* + meson. Our 
predictions for g,j/^ Vc j — 0.681(0.673] and 5d«+d+ 7 = 
— 0.384[— 0.398] obtained from the linear[HO] parame- 
ters fall within the experimental error bars. Our result 
for the coupling constant ratio I 9d *° d °~i I — 4.64(4.59] 

obtained from the linear [HO] parameters is quite com- 
parable with other theoretical model predictions such as 
those 6.32±2.97 [1 and 3.05±0.63 Q from the QCD sum 
rules, 5.54 ± 3.00 9] from the heavy quark effective the- 
ory(HQET), and 4.49 ± 0.96 [HI from the broken-SU(4) 
symmetry by Ml transition. Incidentally, our result for 
the coupling constant ratio I 9b *° b °i I = 0.57(0.57] ob- 

tained from the linear [HO] parameters is the same as 
that from the other QM predictions [2, H, 0] • This result 
is also comparable with 0.64 ±0.51 @ and 0.49 ±0.38 Q 
from the QCD sum rules, and 0.59 ± 0.48 [[| from the 
HQET. 

We show in Fig. [3] our results of decay form factors 
F\/p(q 2 ) obtained from the linear parameters. Since 
the results from the HO parameters are not much dif- 
ferent from those of linear ones, we omit them for 
simplicity. The left (right) panel shows the results of 
charmed(bottomed) vector meson radiative V — > P7* de- 
cays. The solid, dotted, dashed, and dot-dashed lines in 
the left (right) panel represent the form factors for D* + — » 
P+ 7 *(P*+ -> 5+7*), D*° -> D°j*(B*° -> B° 7 *), 
D*+ -> Dfj*(Bf -> P s ° 7 *), and J/i\> -> ^ c7 *(T -> 
?7fc7*) decays, respectively. The arrows in the figure 
represent the zero recoil points of the final state pseu- 
doscalar meson, i.e. q 2 = g^ ax = (My — Mp) 2 . We 
have performed the analytical continuation of the decay 
form factors Pyp(g 2 ) from the spacelike region(g 2 < 0) 
to the physical timelike region < q 2 < q 2 ^^. The cou- 
pling constant gvp-y at q 2 — corresponds to a final 
state pseudoscalar meson recoiling with maximum three- 
momentum |Pp| = {My — Mp)/2Mv in the rest frame 
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TABLE III: Bottomed meson decay constarits(in unit of MeV) obtained from the linear[HO] parameters. 





fs 




Ib s 


/•BS 


fvb 


h 


Linear [HO] 


189[180] 


204[193] 


234 [237] 


250 [254] 


507[897] 


529[983] 


Lattice [24] 


179 ± 18l|* 


196 ± 24±if 


204 ± 16ig° 


229 ± 20t 


- 


- 


QCD [34] 


216 ± 22 


— 


259 ± 32 


— 


— 


— 


\60\ 


1 OQ _|_ 01 




ZoU It OU 








Sum-rules [36] 


210 ± 19 




244 ± 21 








[26] 


203 ± 23 




236 ± 30 








BS [27] 


196 ± 29 


238 ± 18 


216 ±32 


272 ± 20 




498 ± 20 


QM [28] 


155 ± 15 




210 ± 20 








RQM [29] 


189 


219 


218 


251 






Exp. 


229t£t£ [37] 










715 ±5 [19] 




TABLE IV: Coupling constants gvp~/[GeV x ] for radiative 
V — > P7 decays obtained from the linear [HO] parameters. 



Coupling 


This work 


\3] 


[2] 


yj 


Exp. [19] 


9J/ipr} c j 


0.681[0.673] 






0.69 


0.57 ±0.11 


9d*+d+i 


-0.384[-0.398] 


-0.30 


-0.37 


-0.35 


-(0.50 ±0.12) 


9D'°D°-t 


1.783[1.826] 


1.85 


1.94 


1.78 






4.64[4.59] 


6.17 


5.24 


5.08 






-0.167[-0.161] 




-0.17 


-0.13 




9b*+b+i 


1 .311(1.313] 


1.40 


1.50 


1.37 




9b*°b°~, 


-0.749[-0.750] 


-0.80 


-0.85 


-0.78 




1 »B««B» 7 1 
9b' + B + 7 ' 


0.57[0.57] 


0.57 


0.57 


0.57 




9Bf'B0j 


-0.553[-0.536] 




-0.62 


-0.55 




9~fvb~t 


-0.124[-0.119] 






-0.13 





of vector meson. Due to the small kinematic region < 
Q 2 < 9max f° r the bottomed and bottomonium meson de- 



cays, the recoil effects of the final state mesons are quite 
negligible, i.e. F VP (q 2 aax )/g V p 1 w 1. Likewise, we find 

that i ? J/-0r )c (9max)/S'J/V"7o7 K ^D*D B (((max)/ S-DjD s 7 « 1 

for J/ip — » ?y c 7* and Z?* + — ► Dfj* decays. On 
the other hand, we obtain F D ,+ D +(q 2 lax )/g D ,+ D + 1 = 
0.420/0.384 « 1.1 and F J3 .o J3 o(g2 Iut )/ to .o IJ a 7 = 
1.859/1.783 w 1.04 for £>*+ -> L>*+ 7 * and L>*° -> 
D*°7*decays, respectively. The recoil effect, i.e. the dif- 
ference between the zero(g 2 lax ) and the maximum^ 2 = 
0) points, may not be negligible especially for the D* + — > 
D + -f* decay. Fig. [3] also shows the restoration of 
SU(3) flavor symmetry, F D ,+ D ,+ (q 2 )/F D ,+ D ,+ (q 2 ) -> 
1 between charmed and charmed-strange mesons and 
Fb'0B'°{i 2 ) I ' F B *o B *o(q 2 ) — ► 1 between bottomed and 
bottomed-strange mesons in the intermediate and deep 
spacelikc(g 2 < 0) region, where the light quark current 
contribution becomes negligible. 

For a more direct comparison with the available ex- 
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TABLE V: Decay widths and branching ratios for radiative V —* P7 decays obtained from our linear [HO] model parameters. 
We used M Vb = 9353 ± 50 MeV for T -> 1^7 decay. 



Decay mode 


1 [Kevj 




n r r-i Ql 


J/ip — * rj c r ) 


1.69 ± 0.05(1.65 ±0.05] 


(1.80 ±0.10)[1.76± 0.10]% 


(1.3 ±0.4)% 


D *+ _^ D+1 


0.90 ± 0.02(0.96 ±0.02] 


(0.93 ±0.31)[1.00± 0.34]% 


(1.6 ±0.4)% 


P* u -> _D U 7 


20.0 ±0.3(21.0 ±0.3] 




(38.1 ±2.9)% 


-. D s + 7 


0.18 ±0.01(0.17 ±0.01] 




(94.2 ± 0.7)% 


B*+ -> B+7 


0.40 ± 0.03(0.40 ± 0.03] 








0.13 ±0.01(0.13 ±0.01] 






b: u - Bh 


0.068 ±0.017(0.064 ±0.016] 






T — > ^7 


0.045i u ;^[0.042tr Q f 6 ] 


(8.4i£ b )[7.7i^] x 10- 4 





perimental data, we finally calculate the partial decay 
widths from Eq. ([T5]) . In Table [V] we present our 
results for the decay widths and branching ratios to- 
gether with the available experimental data. The er- 
rors in our results for the decay widths and branching 
ratios come from the uncertainties of the experimental 
mass values and experimental mass values plus the full 
widths, respectively. Our results of the branching ra- 
tios Bt(J/4> -> 77 C 7) = 1.80 ± 0.10(1.76 ± 0.10]% and 
Br(D* -► P>+ 7 ) = 0.93±0.31[1.00±0.34]% obtained from 
the linear [HO] parameters are in agreement with the ex- 
perimental data PI, Br(J/ip -> r? c 7) oxp = (1.3 ± 0.4)% 
and Br(P>* -> D+7) exp = (1.6 ± 0.4)% within the error 
bars. For the neutral charmed meson decay, our predic- 
tion T(D*° -> Z?°7) = 20.0 ±0.3[21.0± 0.3] keV obtained 
from the linear [HO] parameters is to be compared with 
other theoretical model results such as 21.69 keV from 
the RQM [|, 14.40 keV from the QCD sum rules Q and 
27. 0± 1.8 keV from broken-SU(4) symmetry by Ml tran- 
sition [39]. For the charmed-strange meson decay, our 
prediction T(D*+ -> D+7) = 0.18 ± 0.01(0.17 ± 0.01] 
keV obtained from the linear [HO] parameters is compa- 
rable with other theoretical model results such as 0.19 
keV from the RQM @ and 0.3 keV and (0.24 ± 0.24) 
keV [| from the HQET. Since the D*° lifetime has not 
been measured yet, we also try to estimate the full width 
for D*° meson using the relation 

Br(D*+ -> £+7) _ T(D*+ -> D+j) r tot (£>*°) 
Br(D*° -» P>° 7 ) T(D*° -> P>° 7 ) r to t(P>* + ) ' 

where we use our predicted decay width T(D*° — > D°-f) 
to extract the full width for D*° . Similarly, we can 
estimate the full width for D* + meson using the same 
method as in the case of D*° meson. Our averaged val- 
ues of the full widths for D*° and D* + mesons obtained 
from the two parameter sets are 

r t ot(D*°) = (55±6)keV, 

r t otOD: + ) = (0.19 ± 0.01) keV, (21) 

respectively, while experimentally only upper limits were 
reported as r(D*°) cxp < 2.1 MeV and r(P>;! + ) CX p < 
1.9 MeV. Some other theoretical model predictions of 
the full widths for D*° and D* + mesons were also re- 
ported as T tot (D*°) = 65.09 keV from the RQM [3] 



and T tot (D*°) = (36.7 ± 9.7) keV and T tot {D*+) = 
(0.24 ± 0.24) keV from the HQET d- 

For B* and B* radiative decays, our results for the 
decay widths T(B*+ -> P+7) = 0.40 ± 0.03[0.40 ± 0.03], 
T{B*° -> P° 7 ) = 0.13 ± 0.01[0.13 ± 0.01], and T{Bf -> 
B°7) = 0.068 ± 0.017[0.064 ± 0.016] obtained from the 
linear [HO] parameters are quite comparable with other 
theoretical model predictions such as T{B* + — > B + j) = 
0.429 keV and T(B*° -> B°j) = 0.142 keV from the 
RQM d, r ( B * + -> B+ l) = (O- 22 ± 0.09) keV and 
T(B*° -> B°7) = (0.075±0.027) keV from the HQET d, 
and T{B*+ -► P+7) = 0.14 keV and T(B*° -► P° 7 ) = 
0.09 keV from the chiral perturbation theory [ll|. Fi- 
nally, for the T — ► 77^7 process, our predictions for the 
decay width and branching ratio obtained from the lin- 
ear [HO] parameters are T(T -> rj b j) = 45±||[42±||] eV 
and Br(T -» r) b -y) = (8.4^ 8 2 6 )[7.7±,i 7 6 ] x 10" 4 , where the 
upper, central, and lower values correspond to Am = 60 
MeV, 110 MeV, and 160 MeV, respectively. For this bot- 
tomonium radiative decay, the decay width T(T — * ?yyy) 
is found to be very sensitive to Am because it is pro- 
portional to (Am) 3 . Our result is to be compared with 
other model predictions such as T(T — > 77^7) = (3.6±2.9) 
eV |l2[ from the nonrelativistic effective field theory 
model, (33. 2 ±0.1) eV [H and 5.8 eV [3 from the RQM. 



V. SUMMARY AND DISCUSSION 

In this work, we investigated the weak decay con- 
stants and the magnetic dipole V ~ * P7 decays of 
heavy-flavored mesons such as (D,D*,D s ,D*,rj c ,J/ip) 
and (B,B*,B s ,B*,ri b ,T) using the LFQM constrained 
by the variational principle for the QCD-motivated ef- 
fective Hamiltonian. The momentum dependent form 
factors Fvp(q 2 ) for V — > P7* decays are obtained in 
the q + = frame and then analytically continued to the 
timelike region by changing qj^ to iq± in the form fac- 
tors. The coupling constants gvpjj which are needed 
for the calculations of the decay widths for V — > P7, 
can then be determined in the limit as q 2 — ► 0, i.e. 
gvp-y = Fyp{q 2 — 0). Our model parameters obtained 
from the variational principle uniquely determine the 
above nonperturbative quantities. This approach can es- 
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tablish the extent of applicability of our LFQM to wider 
ranging hadronic phenomena. 

Our predictions of mass spectra and decay constants 
for heavy pseudoscalar and vector mesons are overall in 
good agreement with the available experimental data as 
well as other theoretical model calculations. Our nu- 
merical results of the decay widths for J/ip — > 7^7 and 
D* + — > D + j fall within the experimental error bars. 
We also estimates the unmeasured full widths for D*° 
and D*+ as r tot (L>*°) = (55 ± 6) keV and T tot (D*+) = 
(0.19 ± 0.01) keV, respectively. Our predictions for the 
branching ratios for the bottomed and bottomed-strangc 
mesons are quite comparible with other theoretical model 
predictions. For the radiative decay of the bottomonium, 
we find that the decay widths T(T — > Vbl) is very sen- 
sitive to the value of Am = My — M Vb . This sensi- 
tivity for the bottomonium radiative decay may help to 
determine the mass of rjt, experimentally. In going be- 
yond the static result to see the momentum dependence 
of the form factor for V — + -P7*, we find that most 
results in the heavy flavored sector stand almost unaf- 
fected from the recoil effects. However, the form factor 
F D ,+ D +(q 2 ) seems to give a non-negligible recoil effect 



about 10% between zero and maximum recoil points, i.e. 
F D , +D+ (g I 2 nax )/ i g Z3 » +D+7 ps 1.1. 

Since the form fator Fvp(q 2 ) of vector meson radia- 
tive decay V — ► Pj* presented in this work is precisely 
analogous to the vector current form factor g(q 2 ) in weak 
decay of ground state pseudoscalar meson to ground state 
vector meson, the ability of our model to describe such 
decay is therefore relevant to the reliability of the model 
for the weak decay. Consideration on such exclusive weak 
decays in our LFQM is underway. Although our previ- 
ous LFQM [HG3 and this analyses did not include the 
heavy mesons comprising both c and b quarks such as B c 
and B* , the extension of our LFQM to these mesons will 
be explored in our future communication. 
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